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Abstract 

The masses of experimentally known highly excited baryons of negative parity supposed to belong 
to the [70, £~] multiplets (I = 1,2,3) of the N = 3 band are calculated in the 1/N C expansion method 
to order 1/N C by using a procedure which allows to considerably reduce the number of linearly 
independent operators entering the mass formula. The numerical fits to present data show that 
the coefficients encoding the QCD dynamics have large, comparable values, for the flavor and spin 
operators. It implies that these operators contribute dominantly to the flavor-spin SU(6) symmetry 
breaking, like for the [70, 1 _ ] multiplet of the N = 1 band. 
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I. INTRODUCTION 



Ml 



The 1/N C expansion method, where N c is the number of colors, |l|-|4| is a powerful and 
systematic tool to study ground state baryons 5-^]. The method is based on the observation 
that, for Nf flavors, the ground state baryons display an exact contracted SU(2iVj) symmetry 
when N c — > oo. At large, but finite N c , this symmetry is broken by contributions of order 
of l/N c , leading to mass splittings. 

Subsequently, efforts have been made to extend this method to excited states. These 
states can be grouped into the so-called excitation bands N = 1,2,3, etc. following a har- 
monic oscillator notation. In this way, one can organize the states into SU(6) x 0(3) 
multiplets. So far the resonances corresponding to the N — 1 band have drawn a particular 
attention, being well known experimentally. It turned out that the problem is more com- 
plicated technically than for the ground states, because these states belong to the SU(6) 
x 0(3) [70, 1~] multiplet, thus have mixed symmetric orbital and flavor-spin parts of the 
total wave function. In such a case the SU(2iVj) symmetry is broken at order 0(N®). The 
standard analysis is based on the se par ation of the system into a ground state core + an 



excited quark, either for Nt 



!6| or for Nf = 3 [17]. A simpler method, avoiding 



ihis separation has been proposed in Ref. [18| for Nf = 2 and extended to Nf = 3 in Refs. 



20]. 



The N = 2 band contains five SU(6) x 0(3) multiplets from which four have a physical 
relevance. The [56', + ] and [56, 2 + ], having symmetric orbital and spin-flavor states, have 
been analysed in Refs. [2l| and [22| respectively, in close analogy to ground states. The 
masses of the multiplet of mixed orbital symmetry [70, £ + ], with £ = and 2, have been 
calculated by extending the ground state + excited quark method to an excited symmetric 



core + excited quark 



23] for Nf = 2. The method has been extended to Nf = 3 in Ref. 



24]. 



The N = 3 band contains eight SU(6) x 0(3) multiplets. In the notation of Ref. [25 1 
these are [56, l - ], [56, 3"], [70', 1~], [70", r], [70,2"], [70, 3"], [20, 1"] and [20, 3~], where 
[70', l - ] and [70", l - ] correspond to radial excitations. This classification provides 45 non- 
strange states ( 1 state N9/2~ , 1 state A9/2~, 5 states N7/2~, 2 states A7/2 - , 8 states 
N5/2~, 4 states A5/2~, 9 states N3/2~, 5 states A3/2", 7 states A^l/2~ and 3 states 
Al/2~). On the other hand in the 1900 MeV - 2400 MeV region only about ten non-strange 
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resonances have been observed so far. The interest in the N = 3 band has been largely 
hindered from a theoretical analysis in the 1/N C expansion method, because of the scarcity 
of experimental data on the one hand, and because of its complex multiplet structure on 
the other hand. To our knowledge the only 1/N C expansion study existing so far is that of 
Ref. [26[ in conjunction with Regge trajectories. The analysis included a series of multiplets 
belonging to the N = 1,2, 3, 4, 5 and 6 bands, in particular the [70, 3~] multiplet of iV = 3. 
The mass operator was reduced to a few terms containing simplified operators, considered 
to capture the main features of the spectrum, but the only term of interest was the leading 
spin-flavor singlet term proportional to N c . 

The N = 4 band has 17 SU(6) x 0(3) multiplets [22] from which only the [56,4+] has 
been analysed in the 1/N C expansion method, being the lowest one and also technically 
simple, due to its symmetric orbital and spin-flavor parts 



The above studies revealed a systematic dependence of the contribution of the dominant 
terms in the mass formula, with the excitation energy, or alternatively with the band number, 
as presented in Ref. 23[. It turns out that the coefficient c\ of the leading spin-flavor 
singlet term, proportional to N c , is raising linearly with N. It was also found that the 
coefficient of the spin-orbit operator having matrix elements of order 0(N®), decreases with 
N and tends to vanish at large excitation energy. The coefficient of the spin-spin term, with 
matrix elements of order 0(N~ ), which brings the largest contribution to the splitting, 
decreases with the excitation energy. Results for the N = 3 band were however absent in 
that analysis. Note that the energy dependence of the mass formula obtained in the 1/N C 
expansion method is remarkably compatible with the energy dependence obtained within 
the framework of quark models with a chromomagnetic hyperfine interaction 29 1. 

Here we present the first systematic attempt towards studying the N = 3 band in the 
1/N C expansion method. We include all mixed symmetric multiplets [70, i~\ (£ = 1,2 and 
3) of the band, therefore more experimental data to analyze. An ultimate aim is to see 
whether or not the results are compatible with the systematic analysis of Ref. [23] and to 
clarify the role of the isospin operator, found so important in the N = 1 band. 

At this stage it is useful to mention that both the symmetric core + excited quark 
procedure and our way of handling the problem of mixed symmetric states are algebraic 
methods in the spirit of the Gell-Mann-Okubo mass formula. There is no radial dependence 
in the picture. The symmetric core + excited quark was originally proposed 10[ as an 
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extension of the ground state treatment to excited states and was inspired by the Hartree 
picture. In this way, in the flavor-spin space, the problem was reduced to the knowledge of 
matrix elements of the SU(2Ny) generators between symmetric states, already known from 
the ground state studies. Accordingly, the wave function was approximately given by the 
coupling of an excited quark to a ground state core of N c — 1 quarks, without performing 
antisymmetrisation. In our approach, all identical quarks are treated on the same footing and 
we have an exact wave function in the orbital-flavor-spin space. As a price, the knowledge of 
the matrix elements of the SU(2Ny) generators between mixed symmetric states is required. 



We have ca 
in Refs. 191 



culated and provided all matrix elements for SU(4) in Ref. 18[ and for SU(6) 



20], by considering an extension of the Wigner-Eckart theorem. While in the 



symmetric core + excited quark procedure the number of terms entering the mass formula 
is excessively large our mass formula has less terms and is physically more transparent. The 
calculated spectrum pointed out the important role of the isospin operator (indirectly of the 
flavor-spin operator, as a part of the SU(2N/) Casimir operator), systematically neglected 
in the symmetric core + an excited quark procedure. We consider that a basic aim of the 
1/N C expansion is to find the most dominant terms with a physical meaning. 

Later on, the symmetric core + an excited quark approach was strongly supported by the 
authors of Ref. [30] . Starting from a general large N c constituent quark model Hamiltonian 
and an exact wave function, they used transformation properties of states and interactions 
under the permutation group Sn c and arrived at an expectation value of the Hamiltonian in 
terms of matrix elements of the approximate wave function of Ref. 13| and operators acting 
either on the symmetric core or on the excited quark. The difference with the algebraic 
methods is that the space degree of freedom enter the discussion. As a consequence extra 
terms appear in the mass formula. They were gathered together to match various quark 
models Hamiltonians. This matching implies constraints on the ratios of different coefficients 
(expressed in terms of radial integrals). The numerical application to the Isgur-Karl model 
was partially successful [3Jj by reproducing the spin-spin term but the tensor term could 
not be reproduced. More work and a deeper understanding is therefore desirable. 

Before describing the 1/N C expansion method let us recall some elements of the history 
of the N = 3 band within the framework of the constituent quark model. An important 
wave of interest has been trigerred by the need of finding an assignment of the -D3s(1930) 
resonance announced by Cutkosky et al. [32]. The quark model calculations of Cutkosky 
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and Hendrick 33j and later on of Capstick and Isgur 34|, incorporating a linear confinement 
and relativistic effects, predicted a mass of about 200 MeV above the experimental value. 
An earlier analysis based on sum rules derived in a harmonic oscillator basis by Dalitz et 



al. 



35] provided a mass of 2088 ± 25 MeV describing this resonance as a a pure [56, f ], 



J = 5/2 state of the N = 3 band, following the suggestion of Cutkosky et al. [32] . A similar 
mass range has been obtained in Ref. 25] in a semirelativistic constituent quark model with 
a linear confinement and a chromomagnetic interaction. The spin independent part of the 
model used in Ref. 25[, which has a linear confinement, makes the [56, 1~] multiplet the 



lowest one among those compatible with the quantum numbers of the 1^35(1930) resonance. 
However this resonance remains an open problem in quark models, inasmuch as its mass is 
about 200 MeV too high above the experimental value. 

As already mentioned, here we are concerned with resonances which can be interpreted as 
members of the mixed symmetric multiplets of the N = 3 band. An important incentive to 
this work was that there is new experimental interest, for example, in the photo-production 
of 1] mesons off protons which suggest a new resonance iV(2070).Di5 which can belong to 
]. Moreover a recent multichannel partial wave analysis including high 



the N = 3 band 



lying resonances, in the so-called fourth resonance region 
a high-lying spin quartet 



37 



suggests the existence of 



JV(2150)3/2-, iV(2060)5/2-, iV(2190)7/2 _ , ^(2250)9/2" 



(1) 



with L = 3. In the following we shall compare this suggestion with our predictions. 

The paper is organized as follows. In the next section we introduce the mass operator 
defined within the 1/N C expansion method. In Section III we present the results of four 
distinct fits of the dynamical coefficients in the mass formula and calculate the mass of 
the fitted resonances obtained from one of these numerical fits. We discuss our results in 
a general context by analogy to results obtained for the N = 1 band and compare our 



interpretation of resonances with that of Refs. 
last section. 



36^38 



38|. Some conclusions are drawn in the 
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TABLE I: Operators and their coefficients in the mass formula obtained from numerical fits. The 
values of Cj and d% are indicated under the heading Fit n (n = 1, 2, 3, 4). 



Operator 




T7i"i_ "1 / TV T 1 r\ 

bit 1 (Mev) 


iit 2 (Mev) 


iit o (Mev) 


tit 4 (Mev) 


Oi = N c l 




ci = 672 ± 8 


Ci = 673 ± 7 


Ci = 672 ± 8 


ci = 673 ± 7 


2 = P-s 1 




c 2 = 18 ± 19 


c 2 = 17 ±18 


c 2 = 19 ± 9 


c 2 = 20 ± 9 


o 3 = —S^ 1 

N c 


c 3 = 121 ± 59 


c 3 = 115 ±46 


c 3 = 120 ± 58 


c 3 = 112 ±42 


1 

O4 = — — 

N c 


T a T a - -^N C (N C + 6) 


c 4 = 202 ± 41 


c 4 = 200 ± 40 


c 4 = 205 ± 27 


c 4 = 205 ± 27 


q = ^_ L i T aQ la 

N c 


cs = 1 ± 13 


c 5 = 2 ± 12 






O = j£ L (2)ij G m G ja 

N c 


c 6 = 1 ± 6 




c 6 = 1 ± 5 




B 1 = -S 




di = 108 ± 93 


d x = 108 ± 92 


di = 109 ± 93 


di = 108 ± 92 


Xdof 




1.23 


0.93 


0.93 


0.75 



II. THE MASS OPERATOR 

When hyperons are included in the analysis, the SU(3) symmetry must be broken and 
the mass operator takes the following general form 39 j 



M = Y J c l O i + Y J d i B i . 



(2) 



The formula contains two types of operators. The first type are the operators Oj, which are 
invariant under SU(Nf) and are defined as 



Oi 



- n (k) n (k) 
-— jU £ ■ u SF , 

c 

where Of^ is a A;-rank tensor in SO(3) and Ocfp a /c-rank tensor in SU(2)-spin. Thus Oi are 
rotational invariant. For the ground state one has k = 0. The excited states also require 
k — 1 and k = 2 terms. The rank k = 2 tensor operator of SO (3) is 



(3) 



L {2)i 3 = Ul\V}- h -jL-L, 



(4) 



2 <■ ' J 3 

which we choose to act on the orbital wave function \£mt) of the whole system of N c quarks 



(see Ref. 



23 1 for the normalization of Z^ 2 ** 4 - 7 ). The second type are the operators Bi which are 
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SU(3) breaking and are defined to have zero expectation values for non-strange baryons. Due 
to the scarcity of data on the hyperons here we consider only one hyperon and accordingly 
include only one of these operators, namely B\ = —S where S is the strangeness. 

The values of the coefficients and di which encode the QCD dynamics are determined 
from numerical fits to data. Table [I] gives the list of Oj and Bi operators together with 
their coefficients, which we believe to be the most relevant for the present study. The choice 
is based on our previous experience with the N = 1 band 20j. In this table the first 
nontrivial operator is the spin-orbit operator O2. In the spirit of the Hartree picture 
generally adopted for the description of baryons, we identify the spin-orbit operator with 
the single-particle operator 

£.« = £*(<).*(*), (5) 

i=l 

the matrix elements of which are of order N®. For simpl icity we ignore the two-body part of 
the spin-orbit operator, denoted by 1/N C (£ ■ S c ) in Ref. 13[, as being of a lower order (there 
the lower case operators £(i) act on the excited quark and S c is the core spin operator). The 
analytic expression of the matrix elements of O2 is given in the Appendix. 

The spin operator O3 and the flavor operator O4 are two-body and linearly independent. 
The expectation values of O3 are simply equal to —S(S + 1) where S is the spin of the 
whole system. They are given in Table [III 

Note that the definition of the operator O4, indicated in Table HI is such as to recover the 

4) by subtracting the quantity (N c + 6)/12. 



ljj for the matrix elements of l/N c (T a T a ) 



matrix elements of the usual l/N c (T a T a ) in SU 
This is understood by using Eq. (30) of Ref. 
extended to SU(6). Then, as one can see from Table HT1 it turns out that the expectation 
values of O4 are positive for octets and decuplets and of order N~ , as in SU(4), and negative 
and of order N® for flavor singlets (see the Appendix for details). 

The operators O5 and 0§ are also two-body, which means that they carry a factor 1/N C 
in the definition. However, as G m sums coherently, it introduces an extra factor N c and 
makes the matrix elements of O5 and 0$ of order N®, as seen from Table HT1 These matrix 
elements are obtained from the formulas (B2) and (B4) of Ref. [3] where the multiplet 
[70, 1~] has been discussed. Interestingly, when N c = 3, the contribution of O5 cancels out 
for flavor singlets, like for I = 1 20]. This property follows from the analytic form of the 



isoscalar factors given in Ref. 



2Q|. 



7 



TABLE II: Diagonal matrix elements of the operators Oi for the [70, i ] multiplets ( I = 1, 2, 3) 
of the N = 3 band. 







Oi 


o. 


o 3 


Oa 


o 5 


O e 


4 JV[70,3" 


h/2 


iV c 


3 
2 


15 

TnZ 


3 


9(iV c + 3) 
4AT C 


75(7V C - 1) 
8Af c 


2 JV[70,3" 


h/2 


N c 


2iV c -3 
2iV c 


3 


3 


9 

2iVc 





4 JV[70,3" 


h/2 


N c 


7 
6 


15 


3 

1n~ c 


7(iV c + 3) 
ANc 


45(iV c - 1) 
8N C 


2 iV[70,3" 


h/2 


N c 


2(27V C + 3) 
3N C 


3 


3 


6 





4 7V[70,3- 


h/2 


N c 


-2 


15 


3 


3(iV c + 3) 
N c 


45(7V C - 1) 
2N C 


2 AT[70',1" 


'h/2 


N c 


2Nc -3 


3 

Tn~ c 


3 


3 
2iVc 





2 N[70', V 


"]l/2 


N c 


2Af c - 3 
3iV c 


3 


3 


3 

~N~c 





2 A[70,3" 


7/2 


N c 


1 

2 


3 


15 


9(N C + 1) 
W c 





2 A[70,2" 


5/2 


N c 


1 

3 


3 

Tn~ c 


15 


3(N C + 1) 
2AT C 





2 A[70,3" 


7/2 


N c 


3 
2 


3 

4We 


2Af c + 3 
4iV c 


3(N C - 3) 
4Af c 






We remind that the SU(6) generators S\ T a and G ia and the 0(3) generators V of Eq . 



act on the total wave function of the N c system of quarks as proposed in Refs. 18], 




201 ] . The advantage of this procedure over the standard one, where the system is 



separated into a ground state core + an excited quark [13j, is that the number of relevant 
operators needed in the fit is usually smaller than the number of data and it allows a better 
understanding of their role in the mass formula, in particular the role of the isospin operator 
O4 which has always been omitted in the symmetric core + excited quark procedure. We 
should also mention that in our approach the permutation symmetry is exact [ls| . 

Among the operators containing angular momentum components, besides the spin-orbit, 
we have included the operators O5 and Oq, to check whether or not they bring insignificant 
contributions, as it was in the N — 1 band. From Table [I] one can see that their coefficients 
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are indeed negligible either included together as in Fit 1 or separately as in Fit 2 and 3. 



TABLE III: Partial contributions and the total mass (MeV) predicted by the 1/N C expansion 
method obtained from Fit 4. The last two columns indicate the empirically known masses and the 
resonance name and status (whenever known). 









Part, contrib. (MeV) 




Total (MeV) 


Exp. (MeV) 


Name, status 






ciOi 


C2O2 


C3O3 


C4O4 










4 iV[70,3" 


h/2 


2018 


29 


140 


51 





2238 ± 46 


2275 ± 75 


Gi 9 (2250)**** 


2 iV[70,3~ 


h/2 


2018 


10 


28 


51 





2107 ± 17 


2150 ± 50 


Gi 7 (2190)**** 


4 iV[70,3" 


Is/2 


2018 


-23 


140 


51 





2186 ±41 


2180 ± 80 


Di 5 (2200)** 


2 iV[70,3~ 


h/2 


2018 


-39 


28 


51 





2058 ± 14 


2060 ± 15 


£>i 5 (2060) 


4 iV[70,3- 


h/2 


2018 


-39 


140 


51 





2170 ± 42 


2150 ± 60 


D 13 (2150) 


2 AT[70',r 


~h/2 


2018 


3 


28 


51 





2101 ± 14 


2081 ± 20 


Di 3 (2080)* 


2 N[70',r 


"]l/2 


2018 


-7 


28 


51 





2091 ± 12 


2100 ± 20 


Si 1 (2090)* 


2 A[70,3~ 


h/2 


2018 


-10 


28 


256 





2292 ± 25 


2200 ± 80 


G 3 7(2220)* 


2 A[70,2" 


h/2 


2018 


-7 


28 


256 





2295 ± 25 


2305 ± 26 


D 35 (2350)* 


2 A[70,3 _ ; 


W/2 


2018 


29 


28 


-153 


108 


2030 ± 82 


2030 ± 82 


G 7(2100)**** 



III. RESULTS AND DISCUSSION 

We have performed four distinct numerical fits of the mass formula ([2]) to the experimental 
data. The corresponding dynamical coefficients c, and di together with the values of Xdof are 



9 



listed in Table HI Fit 1 is made with all operators. Fit 2 and Fit 3 are made by removing one 
operator and Fit 4 is made only the first four operators. It turns out that the contributions 
of angular dependent operators 5 and 6 are negligible but that of the spin-orbit operator, 
which is quite small, remains important. The values of Xdof are g°°d and the error bars of 
the coefficients suggest that the choice of the operators we have made provides a reliable fit. 

As already mentioned, in the 1.9 — 2.4 GeV region the experimental data are rather 
scarce. For this reason, besides the eight resonances provided by the Particle Data Group 



37 



. They do not have 



(PDG) 4(3] we have also included two more, proposed in Refs. 
a status yet. Therefore, in all, we have included ten resonances in the fit from which one is 
a hyperon. 



The experimental masses of four star resonances are from the Summary Table of PDG. 
For the two star resonance -Dis(2200) we took the mass indicated in the Baryon Particle 
Listings of PDG as due to Cutkosky et al. [4l| and for the one star resonance _Di 3 (2080) 
the mass due to Hoehler et al. [421 . The experimental mass of S , n(2090) was taken as the 
average of the masses obtained by Hoehler et al. and Cutkosky et al. in the partial wave 
analysis of the irN scattering. 

The baryon masses obtained from the mass formula ([2]) with the coefficients from Fit 
4 and matrix elements from Table [Til are presented in Table IHH Partial contributions of 
different operators to the total mass are also indicated. One can see that the spin operator 
O3 brings a dominant contribution to the splitting in 4 iV[70, 3~] states and the isospin 
operator is as important, or even more, in the A and the A resonances. In the latter, the 
negative sign of O4 matrix elements helps in lowering the calculated mass close down to the 
experimental value. 

The first observation regarding the multiplet structure is related to the the Gi7(2190) 
four star resonance. To obtain a good fit we had to interpret it as a 2 iV[70, 3~] state. If so, 
Table [ 



Refs. 
of Ref. 



II 
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imp lies that it forms a doublet with the newly suggested -Dis(2060) resonance of 



381 ] . If this resonance is the same as the newly suggested iV(2070)-Di5 resonance 



36] we should be in agreement with the latter authors, who proposed a doublet. 



As a matter of fact we have found out that the inclusion of the resonance 5*11(2090), 
interpreted as a radial excitation belonging to the [70', 1~] multiplet, improves the fit, which 
gives confidence in this interpretation. This resonance appears therefore as the spin-orbit 
partner of L> 13 (2080). 
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Regarding the A resonances, our analysis shows that we have interpreted 1)35(2350) as 
a member of a [70, 2~] multiplet. This is inspired by the quark model results of Ref. 25 1 
where this multiplet, having I = 2, is the highest in the spectrum of the spin- independent 
Hamiltonian with a relativistic kinetic energy and a linear confinement. Such a high value is 
expected to lead to a mass as large as that of the mentioned resonance. This choice suggests 



a kind of agreement between quark models and the present fit, well in the spirit of Ref. [29]. 

The value of the coefficient c\ found in our best fit c\ = 673 ± 7 MeV is smaller than 
ci = 731 ± 17 MeV of Ref. 26) but not far from the estimate c% ~ 640 MeV which can 



be extracted from Fig. 1 of Ref. 



23) where a linear dependence of C\ on the excitation 



energy, or alternatively on the band number N, was found. Note that such an energy 



dependence is reproduced by the formula (29) of Ref. 29j where the compatibility between 
the 1/N C expansion method and semi-relativistic quark models with a linear confinement was 
discussed. This compatibility is confirmed by the present results. The value of c 2 = 20 ± 9 
MeV is practically identical with that obtainable from Fig. 1 of Ref. {23]. 

From the comparison of our results with the "new" resonances reported in Refs. 37 
we can make the following comments. We do not support the interpretation of the -Dis(2060) 
resonance as a member of the quartet (pQ), inasmuch as we interpret this resonance as a 
member of the doublet 2 iV[70,3~]. But £>i 3 (2150) is a member of the quartet 4 iV[70,3~] 
together with G 19 (2250) and D 15 (2200). If this interpretation is valid it remains to find the 
J = 7/2 member, not observed yet. 

IV. CONCLUSIONS 

Using the 1/N C expansion method we have analyzed the multiplet structure of high-lying 
negative parity resonances, located in the 1900 MeV - 2400 MeV region, supposed to belong 
to the N = 3 band in the SU(6) x 0(3) classification. Our results are largely consistent 
with the recent experimental analysis of Refs. 36N38) . A possible future observation of a 
7/2 _ resonance would be of great help in understanding the 4 iV[70,3~] multiplet and the 
structure of the N = 3 band in general. 

The simplified method we have used allows us to include a small number of terms in 
the mass formula and easily identify the most dominant operators to order 0(N~ 1 ). As a 
common feature with the SU(4) and SU(6) analysis of the N = 1 band, we found that the 
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isospin operator O4, neglected in the standard core + excited quark approach, contributes to 
the mass of A's with a coefficient C4 with a magnitude comparable to that of the coefficient 
C3 of the spin operator O3 in N* resonances. In addition the role of the operator O4 is crucial 
in describing the flavor singlet four star resonance A(2100)Go7 included in the fit. 

Future discoveries will help to improve our study and confirm or infirm the present mul- 
tiplet interpretation. 



Appendix A 

We remind that the matrix elements of the spin-orbit operator O2 between states with 
spins S and S' are given by 

(£' 3'J'J' 3 ; I'I' z \l ■ s\£SJJ 3 ; II 3 } P=2 = 

( i 1 



(-iy W+1/2 h'j5 4J MSpi5i' 3h ^(2S + 1)(25' + !)£(£ + 1)(2£ + 1) 



S S' J 



EC-i)^^ 1 ' 11 ^)^- 1 ' 11 ^)^ 1 2 2 \. (ai) 

" \ S r S S' 



The quantities c££ l ' 1] (S), are a short hand notation for the isoscalar factors of the permu- 



tation group of N c quarks, denoted by K([f']p'[f"]p"\[f]p) in Ref. |43j (see also Ref. |44|). 
In that notation p, p\ and p 2 represent the position of the Af c -th quark in the spin-flavor, spin 
and flavor parts of the wave function, of partitions [/], [/'] and [/"] respectively. Actually 
c!p / p~ 1 ' 1 \S) are functions of the spin S and the number A^ c of quarks. By definition, the core 

+e rn d quark wave fmiction fl has p ■ 2 This expression is — to Eq - (A7) of 

Ref. [13J. The correspondence in the isoscalar factors denoted there by c pv is 

c[ N r^(S)^c ^ c [ 2 f- hl] (S)^c 0+ ; C ^(^)^c ++ ; cg^S) ^ c__. (A2) 

We also remind that the matrix elements of the isospin operator O4 as defined in Table 
[U requires the knowledge of the expectation value of the SU(3) Casimir operator T a T a . 
Labelling the flavor states by (A, p) this is 

(T a T a ) = h\ 2 + p 2 + Xp + 3X + 3p) (A3) 

For the flavor states of N c quarks we have 
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W8: a = i, „ = *fi, (t«t« H ^>! 

2l0: A — 3, „=^, (r-n = jV ? + W ° + 45 , 
p 2 x 7 12 

21 X n ^-3 /rriCLrriCL\ A c (A c + 6) - 3 (2A C + 3) 



! 1 : A = 0, /i = o , (T a T a ) 



2 12 
The last case has been discussed in Ref. [isj]. Accordingly the expectation value of O4 is 
negative for flavor singlets. 

The analytic expressions of the matrix elements of O5 and Oq can be found in Ref. 20] 
together with the corresponding isoscalar factors of the SU(3) generators which we do not 
reproduce here. 
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